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Abstract 

The concurrence of two alternate qubits in a four-qubit Heisenberg XX chain 
is investigated when a uniform magnetic field B is included. It is found that 
there is no thermal entanglement between alternate qubits if B is close to zero. 
Magnetic field can induce entanglement in a certain range both for the antifer- 
romagnetic and ferromagnetic cases. Near zero temperature, the entanglement 
undergoes two sudden changes with increasing value of the magnetic field B. 
This is due to the changes in the ground state. This novel property may be used 
as quantum entanglement switch. The anisotropy in the system can also induce 
the entanglement between two alternate qubits. 
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The existence of entanglement shows many interesting properties in quantum sys- 
tems. Its nonlocal quantum correlation has become one of the most valuable resources 
in quantum communication [1-3] and quantum computation [4, 5]. Recently, the con- 
cept of thermal entanglement in solids was introduced and studied in one-dimensional 
anisotropic Heisenberg model [6] . Thermal entanglement in two-qubit Heisenberg XX 
chain was investigated with and without an external magnetic field [7-9]. The entan- 
glement between qubits of the next and next-next neighbors in an open spin chain and 
in multi-qubit Heisenberg model was presented [10-13]. 

In most of the previous investigations, the concurrence of two nearest-neighbor 
qubits is calculated as a measure of entanglement. For a pair of two qubits, the con- 
currence is given by [14,15] 

C = max{Ai - A2 - A3 - A4, 0} (1) 

where the quantities = 1,2,3,4) are the square roots of the eigenvalues of the 
operator 

^ = Pl2K®a|)pt2K®a|) (2) 

in descending order. The values of the concurrence are ranged from zero to one when 
quantum states are changed for unentangled to maximally entangled states. 

The state of the system at thermal equilibrium is represented by the density operator 

p(T) = iexp(-^) (3) 

where Z — Tr[exp{—H / kT)] is the partition function, k is the Boltzmann's constant 
and T is the temperature. Since p(T) represents a thermal state, the entanglement in 
the state is called thermal entanglement [16-17]. 

In this report, a Heisenberg XX model of four-qubit in a linear chain is investi- 
gated when a magnetic field B is included. The pairwise entanglement between alter- 
nate qubit is calculated. The four-qubit XX M Heisenberg model is described by the 
Hamiltonian 

4 4 

HxXM = J E (^n + ^^n+l) + ^ E < (4) 

n=l n=l 
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where cr^ are the raising and lowering operations, B is the external magnetic field and 
perpendicular to the chain, J is the strength of interaction. The value of positive and 
negative J corresponds to the antiferromagnetic and ferromagnetic cases respectively. 

The eigenvalues and eigenstates of the Hamiltonian in Eq. (4) can be calculated 
analytically. In a four-qubit Heisenberg chain with periodic boundary conditions, the 
eigenvalues are given by 



Eq = -4S, Ei = 2J - 2B, 

E3 = -2 J - 2B, E5 = 2V2J, 

E7 — Eg — Eg — Elo — 0, 

£^12 = Eu = 2B, Ei3 — —2 J + 2B, 



E2 — E4 — — 2B 
Ee = -2V2 J 
Eu = 2J + 2S 

£15 = 4S 



(5) 



and the corresponding eigenstates can be explicitly expressed as 
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If the concurrence of two alternate qubits is considered, the reduced density matrix 
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Pi3 = Tr24p{T) can be given by 
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in the basis {|00) , |01) , |10) , |11)}. Where u,v,y a^ie 

+ e(2<^+2«)/3 + e-(2^-2S)/3 + cosh2V2Jp) + e2^^ + e^^/^ 
^ 1 (l + e-(^J+2B)p ^ g(2j-2B)/3 ^ cosh2V2J(3) + + 
y = 1 (-1 + cos/i (2 J + 2S) /3 + cos/i (2 J - 2S) /3 + cosh2V2Jp) - cosh2Bp 
The partition function of the system is 

Z = 4 (1 + cosh2Bp)+2 [cosh4:Bp + cosh (2 J + 2B)P + cosh (2 J - 2S) /3 + cosh2V2jp] 



(9) 



with /9 = -j^. In the foUowing calculations k is set to be 1.0. 

From the Eqs. (1), (2) and (7), the concurrence can be obtained 



C = -|max(| 



y\ 



uv, oj 
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Except the concurrence C of two qubits, the global entanglement Q of many-qubit 
pure states also needs to be considered. The global entanglement Q is introduced by 
[18-20] 



1=1 



(11) 



where the ^-concurrence of ICi means the entanglement between the qubit i and the 
other qubits and can be expressed as 



IQ = ^2[1-Tr{pj)] (12) 

The concurrence C as functions of the magnetic field B and the temperature T is 
plotted in Fig. 1. The strength of interaction J is chosen to be 1.0. Fig. 1(a) is a 



three-dimensional plot of the concurrence C as functions of B and T. Prom Fig. 1(a), 
it is clear that there is a two-peak structure in C. The two peaks appear symmetrically 
at two sides of B. There is no entanglement between alternate qubits at i? = 0. The 
entanglement is also independent of T [10]. This can be understood since B — Q 
the ground state will be j'^e), which is unentangled between alternate qubits. When 
the temperature 7 = 0, the ground state energy as a function of the magnetic field 
B is changed from to £^3, and then to Eq for antiferromagnetic case. While the 
ground state energy is changed from E^ to £^11, and then to -E15 for ferromagnetic case. 
Therefore there should be a change in the concurrences C because of the change in 
the ground state. It can be seen that the increase of magnetic field B cannot induce 
entanglement when the temperature T is very high. When the temperature T is very 
low, the entanglement is increased with increasing value of the magnetic field \B\ to a 
maximum value. Then it is decreased and finally disappeared. The maximum value of 
the entanglement decreases when the temperature T is increased. 

The contour map of the concurrence C as functions of the magnetic field B and the 
temperature T is plotted in Fig. 1(b). Four contours of C = 0.5, 0.3, 0.1, are shown 
respectively. Beyond the contour C = 0, the entanglement is equal to zero. It means 
that there exists a critical temperature Tc. From the curve of C = in Fig. 1(b), 
it can be seen that the critical temperature Tc depends on the magnetic field B. If 
B < 0.09, the concurrence C is always zero no matter the temperature T is increased 
or decreased. When 0.09 < B < 0.41 or B > 1.0, there are two critical temperatures 
of Tc. When T is either above the lower part of the curve of C = or below the upper 
part of C = 0, the concurrence C is always greater than zero. When 0.41 < i? < 1.0, 
there is a single value of the critical temperature Tc. When the temperature T is above 
the curve of C = 0, the entanglement is vanished no matter the magnetic field B 
is increased or decreased. For four-qubit nearest-neighbor XXM model, the critical 
temperature is independent of the magnetic field B [11]. 

The thermal entanglement of alternate qubits for J = —1.0 is also studied. The 
result is almost the same as that shown in Fig. 1. This means that the entanglement 
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exits for both antiferromagnetic and ferromagnetic cases. 

The concurrence C, the global entanglement Q, and the ^-concurrence ICi are plot- 
ted in Fig. 2 as a function of the magnetic field B. The concurrence C of different 
qubits is plotted in Figs. 2(a) and 2(b) when the temperature T is varied with T — 0.01, 
0.1, and 0.5. Fig. 2(a) is a plot of C for alternate qubits. From Fig. 2(a), it can be 
seen that the shape of C is like a square wave for low temperature of T = 0.01. The 
concurrence C keeps zero until B is increased to 0.41. Then the concurrence maintains 
a maximal value of C = 0.5 until it drops to zero aX B = 1.0. In the limit of T — > 0, 
one has 

limC = ISI < (V2- 1)|J| 

lim^C=| (V2- 1)|J| < |S| < |J| (13) 

limC = \B\>\J\ 

This can be understood as follows. When \B\ > \J\ and \B\ < {V^ — 1)| J|, the ground 
states are the unentangled state j'^o) and ji/'e) respectively. While for {^/2 — 1)| J| < 
\B\ < \J\, the ground state is the maximally entangled state 1^3) of the two alternate 
qubits. When the temperature T is increased to 0.1, and 0.5, the shape of C is changed 
from square wave to a single peak. The maximal value of the concurrence C is decreased 
when T is increased. For finite temperatures, this definite ground state structure 
is smoothed out by the partition of higher states and therefore the concurrence C 
becomes smaller with higher temperatures. To compare this with that of the two 
nearest-neighbor qubits in the four-qubit XX M model, the concurrence C of the two 
nearest-neighbor qubits is plotted as a function of B in Fig. 2(b) with the same 
condition. For low temperature of T = 0.01, the entanglement keeps a constant value 
of C = 0.46 until it drops to a dip. It seems that the dip is due to the crossing of 
energy level at the point oi B = (\/2 — 1) [11]. Then C maintains a maximal value 
of C = 0.5 until it drops to zero at B = 1.0. When the temperature T is increased 
to 0.1, the shape of C is changed to two peaks with almost the same dip. When T 
is increased to 0.5, the dip disappears. There is only a single peak a,t B — 0. From 
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Figs. 2(a) and 2(b), it can be seen that the effects of temperature on C is much 
stronger for alternate qubits than that for nearest qubits. When T is increased to 0.5, 
the curve of C in Fig. 2(b) is much higher than that in Fig. 2(a). It seems that 
the temperature affects the entanglement between weakly interacting alternate qubits 
stronger than that for strongly interacting nearest qubits. The global entanglement 
Q and the i-concurrence ICi are plotted in Fig. 2(c). Due to the symmetry of the 
eigenstates in Eq. (6), the values of i-concurrence are the same. When < S < 0.41, 
ICi = IC2 = IC3 = IC4 = 1.0. The value of i-concurrence is the same as that of 
Q. The two-qubit concurrence C is C12 — Cu — C23 — C34 = {2^/2 — l)/4 and 
Ci3 — C24 — 0. These are shown in Figs. 2(a) and 2(b). Although the entanglement 
of three and four qubits cannot be discriminated, the additional entanglement of three 
and four qubits can be described by 1 -ECf^ = 0.16 [18, 19]. When 0.41 <B< 1.0, 
ICi = IC2 = ICs = IC4 = a/3/2. The global entanglement Q and the two-qubit 
concurrence C satisfy the relation oi Q — \ Y^ijCfj [18, 19]. There is no additional 
entanglement of three and four qubits. When B > 1.0, there is no entanglement at all. 
All the values of C, Q and ICi equal to zero. It is very interesting to note that in the 
low temperature limit the entanglement between two alternate qubits of the four-qubit 
XXM Heiscnbcrg model undergoes two sudden changes when the magnetic field B 
is increased. This novel property may be used as quantum entanglement switch in 
quantum computing and quantum communications. 

For a more general model of four-qubit Heisenberg XX chain, the anisotropic con- 
tribution needs to be considered. If the anisotropy is included, the Hamiltonian of the 
system can be written as 

J A ^ 

HxxzM = HxxM + — KK+i (14) 

^ n=l 

where A is the anisotropy parameter, Hxxm is given by Eq. (4). The system reduces 
to XX model when A = and the isotropic XXX model when A = 1.0. 

The concurrence C of alternate qubits is plotted in Fig. 3 as functions of the tem- 
perature T, the magnetic field B, and the anisotropic parameter A when the strength 
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of interaction J is 1.0. The concurrence C is plotted as functions of T and A in Fig. 
3(a) when S = 0.5. Prom Fig. 3(a), it is seen that the concurrence C is increased first, 
then reached a maximum value, and finally decreased when A is increased from —0.4 
to 1.0. For A > 0, the concurrence C is monotonically decreased when T is increased. 
The concurrence C is plotted as functions of B and A in Fig. 3(b) when T — 0.2. From 
Fig. 3(b), it is seen that the value of C is kept zero for A > when B — The peak 
in C appears at A = —0.5. It is found that the anisotropy in the Heisenberg model 
can induce the entanglement between alternate qubits even for 5 = 0. If -B > 0, the 
height of the peak in C is increased and the position of the peak is shifted to larger 
values of both B and A. From Fig. 3, it is seen that the curve of C is asymmetric 
about A. 

In conclusion, in this report the entanglement between two alternate qubits of a 
four-qubit Heisenberg XX model is investigated. There is no thermal entanglement 
between alternate qubits of a four-qubit Heisenberg XX model when B is very small. 
However, when the magnetic field \B\ is increased, it can induce entanglement in the 
XX model both for the antiferromagnetic and ferromagnetic cases. The square wave 
like shape appeared in the concurrence C may be used as quantum entanglement 
switch. It is very interesting to find that the temperature affects the entanglement 
much stronger for weakly interacting alternate qubits than that for strongly interacting 
nearest qubits. The anisotropy in the Heisenberg model can also induce entanglement 
between alternate qubits. 

Acknowledgement 

It is a pleasure to thank Yinsheng Ling, Jianxing Fang, and Xiang Hao for their 
many helpful discussions and calculations. 



8 



REFERENCES 

[I] C. H. Bennett, G. Brassard, C. Crepeau, R. Jozsa, A. Peres, and K. Wooters, 
Phys. Rev. Lett. 70, 1895 (1993). 

[2] A. K. Ekert, Phys. Rev. Lett. 67, 661 (1991). 

[3] M. Murao, D. Jonathan, M. B. Plenio, and V. Vedral, Phys. Rev. A59, 

156 (1999). 
[4] B. E. Kane, Nature 393, 133 (1998). 

[5] C. H. Bennett and D. P. Divincenzo, Nature 404, 247 (2000). 

[6] M. C. Arnesen, S. Bose, and V. Vedral, Phys. Rev. Lett. 87, 017901 (2001). 

[7] X. G. Wang, Phys. Rev. A64, 012313 (2001). 

[8] G. L. Kamta and A. F. Starace, Phys. Rev. Lett. 88, 107901 (2002). 

[9] Y. Sun, Y. G. Chen, and H. Chen, Phys. Rev. A68, 044301 (2003). 

[10] L. F. Santos, G. Rigohn and C. O. Escobar, Phys. Rev. A69, 042304 (2004). 

[II] X. G. Wang, H. Fu, and A. I. Solomon, J. Phys. A34, 11307 (2001). 
[12] X. G. Wang, Phys. Rev. A66, 034302 (2002). 

[13] X. Q. Xi, W. X. Chen, S. R. Hao, and R. H. Yue, Phys. Lett. A300, 
567 (2002). 

[14] S. Hill and W. K. Wootters, Phys. Rev. Lett. 78, 5022 (1997). 
[15] W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

[16] M. A. Nielsen, Ph. D. thesis. University of New Mexico, 1998, quant-ph/0011036. 
[17] M. C. Arnesen, S. Bose, and V. Vedral, Phys. Rev. Lett. 87, 017901 (2001). 
[18] J. Endrejat and H. Biittner, quant-ph/0406006. 

[19] U. Glaser, H. Biittner, and H. Fehske, Phys. Rev. A68, 032318 (2003), 

quant-ph/0305108. 
[20] G. Brennen, QIC 3, 619(2003), quant-ph/0305094. 



9 



FIGURE CAPTIONS 
Fig. 1. 

The concurrence C is plotted as functions of magnetic field B and temperature T 
when J — 1.0. 

(a) . The curve of C as functions of B and T. 

(b) . The contour map of C as functions of B and T. 
Fig. 2. 

The concurrence C, the global entanglement Q, and the i-concurrence ICi are plot- 
ted as a function of B when J = 1.0. For the concurrence C of (a) and (b), the curve 
is plotted when T — 0.01, 0.1, and 0.5 (from top to bottom). 

(a) . The curve C of the alternate qubits. 

(b) . The curve C of the nearest- neighbor qubits. 

(c) . The curves of Q and ICi are plotted when T = 0.01. 
: The curve of Q; : The curve of ICi. 

Fig. 3. The concurrence C between alternate qubits is plotted as functions of T, 
B, and A when J = 1.0. 

(a) . The concurrence C as functions of T and A when S = 0.5. 

(b) . The concurrence C as functions of B and A when T = 0.2. 
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Fig. 1 
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